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Abstract 

An accidental degeneracy of resonances gives rise to a double pole in the scattering 
qq | matrix, a double zero in the Jost function and a Jordan chain of length two of general- 

ized Gamow-Jordan eigenfunctions of the radial Schrodinger equation. The generalized 
Gamow-Jordan eigenfunctions are basis elements of an expansion in bound and resonant 
, energy eigenfunctions plus a continuum of scattering wave functions of complex wave 

number. In this biorthonormal basis, any operator f(Hr) which is a regular function of 
the Hamiltonian is represented by a complex matrix wich is diagonal except for a Jordan 
block of rank two. The occurrence of a double pole in the Green's function, as well as 
the non-exponential time evolution of the Gamow-Jordan generalized eigenfunctions are 
C$ ' associated to the Jordan block in the complex energy representation. 

^ { - Keywords: Non-relativistic scattering theory; Multiple resonances; Resonance reac- 

tions; Berry's phase. 
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1 Introduction 

Lately, the interference effects of resonances, the crossing and anticrossing properties of the 
energies and widths of two unbound levels and the occurrence of a double pole of the scattering 
matrix have aroused a great deal of interest. Some interesting examples of interfering unbound 
two level systems are the T = 1,T = 0, J n = 2+ doublet in 8 Be[@, §, [|, the T = 1, T = 
doublet of p and ui messons and the a — K s doublet of neutral sigma and K mesons |2|, [|, [|, || |7|]. 
A variety of widely differing systems where double poles can occur have been identified, such 
as autoionizing states in complex atoms |J and atomic states in intense laser fields||, |10[] . The 
problem of the degeneracy of resonances also arises naturally in connection with the Berry phase 
of resonant states JTT], [l^, [13|, [14j] which was recently measured by the Darmstadt group [pf 
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Some examples of simple quantum mechanical systems with double poles in the scattering 
matrix have been recently described. Vanroose et al.,|l| examined the formation of complex 
double poles of the S-matrix in a two channel model with square well potentials. Recently, 
Hernandez et al.,[17 investigated a one channel model with two spherical concentric cavities 
bounded by 5— function barriers and showed that a double pole of the S"— matrix can be induced 
by tuning the parameters of the model; Vanroose generalized this model to the case of two finite 
width barriers|TB[. The formal theory of multiple pole resonances and resonant states in the 



rigged Hilbert space formulation of quantum mechanics was developed by Bohm et al.,[19] and 
by Antoniou et al.,|ED . 



In the present paper, we deal with the problem of multiple poles of the scattering matrix 
and the generalized complex energy eigenfunctions associated with them in the framework of 
the theory of the analytic properties of the radial wave functions. 

The plan of this paper is as follows. In sections 2 and 3, we introduce some basic concepts 
and fix the notation by way of a short reminder of resonances and resonant states in the theory 
of the analytic properties of the radial wave functions. Sections 4 and 5 are devoted to a short 
discussion of the no-crossing rule for bound states and its non applicability to resonant states. 
In section 6, we show that a double pole of the scattering wave function (double zero of the Jost 
function) is associated to a chain of length two of Gamow- Jordan generalized eigenfunctions 
and derive explicit expressions for this generalized eigenfunctions in terms of the outgoing wave 
Jost solution, the Jost function and its derivatives evaluated at the double pole. We also 
show that the Gamow- Jordan generalized eigenfunctions in the Jordan chain are elements of a 
complete set of states containing the real (bound states) and complex (resonant state) energy 
eigenfunctions plus a continuum of scattering wave functions of complex wave number. In 
section 7 we derive expansion theorems (spectral representations) for operators f(Hr) which 
are regular functions of the radial Hamiltonian H^f 1 and show that, in this basis, the operator 
f(Hr ) is represented by a complex matrix which is diagonal except for a Jordan block of rank 
two associated to the double zero of the Jost function and the corresponding Jordan chain of 
generalized Gamow- Jordan eigenf unction. We give the normalization and orthogonality rules for 
the generalized eigenfunctions in the Jordan chain associated to the double pole of the Green's 
function in section 8. We end our paper with a summary of results and some conclusions in 
section 9. 



2 Regular and physical solutions of the radial equation 

The non-relativistic scattering of a spinless particle by a short ranged potential v (r) is described 
by the solution of a Schrodinger equation. When the potential is rotationally invariant, the 
wave function is expanded in partial waves and one is left with the radial equation 



d 2 (j>e (Ayr) 
dr 2 



+ 



k 



v (r) 



(1) 



As is usually done when discussing the analytic properties of the solutions of (1) as functions of 
k, rather than starting by defining the physical solutions (k,r), we define the regular and 
irregular solutions of (]!]) by boundary conditions which lead to simple properties as functions 
of k. The regular solution (pi(k,r) is uniquely defined by the boundary condition |21 
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lim(2£ + l)!!r-* - V*(fc,r) 



r~*0 



(2) 



0^(A;, r) may be expressed as a linear combination of two independent, irregular solutions of ([]]) 
which behave as outgoing and incoming waves at infinity, 



Mk,r) = Uk- 1 - 1 [h{-k)h{k,r) - (-iy fi (k)M-k,r)] , 

where fe(—k,r) is an outgoing wave at infinity defined by the boundary condition 

lim exp(—ikr)f£(—k,r) = (+i) e 

r— too 

and fg(k,r) is an incoming wave at infinity related to fe(—k,r) by 



(3) 



(4) 



fi(k,r) = (-l) t ft(-k,r) 

for k real and non-vanishing. 

The Jost function fe(—k) = fe(—k, 0) is given by 



(5) 



fl(-k) = (-lYk e W[f e (-k,r),MKr)} 



(6) 



where W[f, g) = fg' — f'g is the Wronskian. The Jost function fe(—k), has zeroes (roots) on 
the imaginary axis and in the lower half of the complex k plane. 

When the first and second absolute moments of the potential exist, and the potential de- 
creases at infinity faster than any exponential (e.g. if v(r) has a gaussian tail or if it vanishes 
identically beyond a finity radius) the functions fi(—k), <^(/c,r), and k e fe(—k,r), for fixed 
r > 0, are entire function of fc|[21j|. 

Therefore, the derivatives of these functions with respect to the wave number k exist and 
are entire functions of k for all finite values of k in the complex fc-plane. 

The differential equations satisfied by the derivatives of the functions 4>e(k,r) and fi(—k,r) 
with respect to k are obtained from ([I]) taking derivatives with respect to k in both sides of 
the equation, 



d 2 (p e (k,r) 
dr 2 



+ 



— v(r) 



4>i(k,r) = -2A;0 £ (fc,r) 



(7) 



d 2 (f) e (k,r) 
dr 2 



+ 



k - 



v(r) 4>t(k, r) = -Ak<j>t(k, r) - 2<j) t (k, r), 



(8) 



in (0) and (|8|) we have used the notation <fie(k,r) = d(pe(k,r)/dk. Similar expressions are valid 
for the derivatives with respect to k of the outgoing wave solutions fi(—k,r). 
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The scattering wave function ipg(k, r) is the solution of equation ([]]) which vanishes at the 
origin and behaves at infinity as the sum of a free incoming spherical wave of unit incoming 
flux plus a free outgoing spherical wave, 



4 + \k,0) = 



(9) 



and 



lim \ ^ +) (k, r) - h>- } (k, r) - St{k)tif> (k, r) 



;(-) 



(+) 



o. 



(10) 



In this expression U~\k,r) and h[ + \k,r) are Ricatti-Hankel functions that describe incoming 
and outgoing waves respectively, Se(k) is the scattering matrix. 

Hence, the scattering wave function ip^ (k, r) and the regular solution are related by 



^ +) (A;,r) 



and the scattering matrix is given by 



3 t (k) 



k e+1 <j) e {k,r) 
fi(-k) 



M-k) ' 



(11) 



(12) 



The complete Green's function for outgoing particles or resolvent of the radial equation 
may also be written in terms of the regular solution <fie(k, r) and the irregular solution fe(—k, r) 
which behaves as an outgoing wave at infinity 



G { +\k;r,r') 



] Wk) 



(13) 



3 Bound and resonant state eigenfunctions 

Bound and resonant state energy eigenfunctions are the solutions of ([!]) which vanish at the 
origin 



u n£ {k n ,0) = 0, 
and at infinity satisfy the boundary condition 



(14) 



lim 



1 du ne (k n ,r) 



r->oo [u ne (k n} r) dr 
where k„ is a zero of the Jost function, 



ik r , 



0. 



(15) 
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M-K) = 0. 



(16) 



From equations ([I]) and (^) we verify that all roots (zeroes) of the Jost function are associated 
to energy eigenfunctions of the Schrodinger equation. 

Bound state eigenfunctions are associated to the zeroes of fe(—k) which lay on the positive 
imaginary axis k 2 s = —k 2 s < 0, while resonant or Gamow state eigenfunctions are associated to 
the zeroes of the Jost function which lay in the fourth quadrant of the complex fc-plane. 

From (0), (01) and (|T6|), bound states and Gamow or resonance eigenfunctions are related 
to the regular solution (pi{k,r) by 



u n e(k n , r) = N n( l 4> e (k n , r) (17) 

where N n £ is a normalization constant. Due to the vanishing of fe(—k n ), (f)g(k n ,r) is now 
proportional to the outgoing wave solution, k n ,r), of (|I]). Hence, 

u ne (k n ,r) = N^-^L-f e (k n )f e (-k n ,r). (18) 

This expression shows, in a very explicit way, that the Gamow state eigenfunctions u n i(k n , r) 
with k n = K n — ij n and n n > 7 n > 0, are solutions of (|l|) which vanish at the origin and 
asymptotically behave as purely outgoing waves which oscillate between envelopes that increase 
exponentially with r, the corresponding energy eigenvalues S n are complex with Re S n > Im£ n . 

The bound state eigenfunctions u s i(k s , r) are also solutions of (|I|) which satisfy the boundary 
conditions (11) and ([15|), but, in this case, k s = in s with, k s > 0, which means that asymptoti- 



cally the outgoing wave of imaginary argument, fe(—k s ,r), decreases exponentially with r and 
the energy eigenvalue £ s is real and negative. 



4 The no-crossing rule for bound states. 

In the case of bound states, the normalization constant is related to the derivative of the Jost 
function evaluated at k s and it may also be expressed as a normalization integral. The zero 
of the Jost function is on the positive imaginary axis, and the bound state eigenfunction is 
cuadratically integrable (for time reversal invariant forces 4>z(iK s ,r) is real). R.G. Newton gives 
the following expression ||21|| 

N 2 si = ^TT) (^d^) k fi{ks) = 1^ mKr)?dT - (19) 

Since the normalization integral is positive and the function fe{k) is regular at k s = in s , 
the derivative of the Jost function evaluated at k s = in s cannot vanish. Therefore, the zero 
of ft{—k) at k s = ik s must be simple. The corresponding pole in G# (k;r, r'), ip^ + \k,r) and 
Si(k) must also be simple. 

It follows that, in the absence of symmetry, the real, negative energy eigenvalues of the 
radial equation for a one channel problem cannot be degenerate. 
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5 Crossing of resonant states 



In the case of a resonant state, the zero of the Jost function fe(—k) lies in the fourth quadrant 
of the complex fc-plane, 



k n — K n ^Ini (20) 

with K n > 7„ > 0. 

The resonant or Gamow eigenfunction <f>i(k n ,r) is an outgoing spherical wave of complex 
wave number k n and angular momentum i. Therefore, for large values of r, (pf (k n ,r) oscillates 
between envelopes that grow exponentially with r. Hence, the integrals over r must be properly 
defined. This may be done by means of a gaussian regulator and a limiting procedure^] T. 
Berggrenp3|, |24| gives the following expression 



1 



2(1+1) 



dfi(-k) 
dk 



fe(k n ) 



k'n 



lim 



exp(— ur 2 )(l)1{k n , r)dr 



(21) 



The integral in the right hand side is a complex number and it may vanish. 

Since fe(k n ) has no zeroes in the lower half of the complex /c-plane, the left hand side of 
equation (|2T| ) vanishes only when (dfe(—k)/dk)k n vanishes. Then, we have two possibilities, 

i) When (dfz(—k)/dk)k n does not vanish, f(—k), has a simple zero at k = k n , the integral 
in the right hand side of equation (|21| ) does not vanish and the normalization constant, N 1 
occurring in 

ii) When 



nt- 



is given by (f2~T|). 



dfi(-k) 
dk 



(22) 



k'n 



the integral in the right hand side of (ET1) vanishes 



lim / exp(— i>r 2 )(f)1{k n , r)dr = (23) 
^° Jo 

and the Jost function fp(—k) has a multiple zero at k = k n . In this case, the Green's function 
G^\k; r, r'), the scattering wave function ipj, (k,r) and the scattering matrix Se(k) have a 
multiple pole at k = k n . The normalization constant of the Gamow eigenfunction is no longer 
given by @. 

Furthermore, it will be shown below that when fe(—k) has a multiple zero ( a multiple 
resonant pole of rank r in G\ (k; r, r'), ip^ + \k, r, ) and Se(k)) the corresponding complex energy 
eigenvalues are degenerate even in the absence of symmetry. That is, the no-crossing rule does 
not hold for resonant eigenstates. 
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6 Completeness and the expansion in complex resonance 
energy eigenfunctions 

In this section, it will be shown that associated to a double zero of the Jost function (double pole 
of the scattering wave function ip^, + \k, r), the Green's function G\ (k; r, r') and the scattering 
matrix Se(k)) there is a chain of generalized Gamow- Jordan eigenfunctions which together 
with the bound state and resonant state eigenfunctions form a biorthonormal set which may 
be completed with a continuum of scattering wave functions of complex wave number. 

Given two square integrable and very well behaved functions $(r) and x{ r ) which decrease 
at infinity faster than any exponential, the completeness of the orthonormal set of bound state 



and scattering solutions of the radial Schrodinger equation [^1|] allows us to write 



7T 

where < > is the standard Dirac bracket notation 



s bound states 

+ -/ <^\4 + \k r )><4 + \k')\x>dk' (24) 



< $|x >= / §*(r)x(r)dr. (25) 
Jo 

We shall assume that the Jost function fe(—k) has a double zero at k = k m in the fourth 
quadrant of the complex A;'— plane, all other zeroes of fe(-k') in that quadrant being simple. 
Then, the scattering function ip^(k',r) as function of A;'— complex, has one double resonance 

pole at k! = k m and simple resonance poles at k = k n , n = 1,2 m — l,m + 1..., all in the 

fourth quadrant of the complex A;'— plane. The function if)( + >* (k' , r) is regular and has no poles 
in the lower half of the A;'— plane. 

In order to make explicit the contribution of the resonant states to the expansion in eigen- 
functions, the integration contour in the second term in the right hand side of ( |2~3] ) is deformed 
as shown in Fig. 1. 

When the deformed contour C crosses over resonant poles, the theorem of the residue gives 



< &\x > 



Y < $\v s/ >< v s/ \x > 

s bound states 

v-^ r2 

+ > 2-iriRes - 

all resonance poles 



<^\4 + \k')x^\k')\ X > 



+ - <$\4 + \k')><4 + \k')\ X >dk'. 
7T Jc 

The residues may be readily computed from equations ([□]) and (|26|). 
When fe(-k') has a simple zero B>t k — , 



(26) 
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/\ Imk 



Rek 




Figure 1: Integration contour C in the complex A;'— plane. 



where 



and, 



2TriRes 

AiRes 

1 



<<S>\4 + \k>) ><^ + >(k')\ X > 



(+)/ 



7T 



<$\Mk') ><Mk')\x> k' 2{e ' +l) 



fijkn) ( df e (-k') 



< $M*0 > 



(k' - k 



dft(-k') 
dk' 



fi(k') 



k 1 — kn 



<Mk')\x> 



<Mk')\x> 



lim 



lim 



$*(r)4> e (k',r)dr, 



4> £ (k',r)x{r)dr. 



(27) 



(2? 



(29) 



since (pi(k',r) is real and bounded for k! real, the integrals in (|28| ) and ( p9|) exist. 

Furthermore, since <f)e(k n ,r) is an outgoing wave which oscillates between envelopes that 
grow exponentially at infinity and $(r) and x{ r ) are very well behaved functions of r that 
decrease at infinity faster than any exponential, the integrals of the products $*(r)0£(fc n , r) 
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and 4>£{k n ,r)x{r) also exist, and we may take the limit indicated in the right hand side of 



equations fl28|) and (|29|) under the integration sign. 
Therefore, 



-<<S>\4 + \k>) ><^>(k')\ X > 

L7T 



2mRes 
where the notation means 



.(+)/ 



=< §\u nl {k n ) >< u nl {k n )\x > (30) 



POO 

<$\u ne (k n )>= §*(r)u nl (k n ,r)dr 
Jo 



(31) 



and 



<u nl {k n )\x>= \ u n e(k n ,r)x(r)dr, 



(32) 



The Gamow eigenfunction or normal mode, u n i(k ni r), is given by (|18l) and the normalization 
constant N n £ is given by 



iAk 



1 f.(k )( df ^- k 'h 



(33) 



in agreement with Berggren's result given in equation (|2T|). 

When the Jost function fe(—k') has a double zero at k' = k m , ip[ + \k',r) has a double pole 

ctij k — kf-n i 



<j) £ {k',r)k'( e+ V 



(k 1 k m )^g£ m [k') 

the function ge m {k') is regular at k! = k m and may be expanded as 



(34) 



9lm{K)- 2 \ ^ ) km + R W Km){ ^ J fcm 4 
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dk' 3 



(35) 



with 



d 2 fe(-k>) 



dk' 



2 



^ 0. (36) 



The function ip£*(k',r') is regular at k' = k m , since fe(k') has no zeroes in the lower half of 
the complex A;'— plane, 
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4 +> (k',r') 



i r i\t,'(e+i) 



fi(k') 



(37) 



Thus, the residue of the term (2/n) < \k') >< 4 \k')\x > at the double pole in k' = k r , 

is obtained from the Cauchy integral formula as 



2mRes 
AiRes 
4z 



-<^\4 + \k') ><^>(k')\ x > 



(+)/ 



7T 



< $\Mk') >< Mk')\x > k' 2{e+1) 
(k 1 - k m ) 2 g im (k>)f e (k>) 

d / < ®\Mk') >< Mk')\x > k' 2(e+1) 
dk\ 9i m {k')fe{k') 



(3* 



After computing the derivative indicated in ( |3~8D and rearranging some terms, we obtain 



2iriRes 



< 



$\4 + \k')><4 +) (k')\ x >} k , 



L7T 



Af 2 



< $\fa(km) >< M k m)\x > + < ®\Mkm) >< Mk m )\x > 



(39) 



where, according to QT7D, 4>£(k m ,r) is the non-normalized Gamow eigenfunction, and (f>i(k m ,r) 
is a generalized Gamow- Jordan eigenfunction or abnormal mode given by 



9l{k m ,r) = — — \-C e (k m )(pe(k m ,r), 



(40) 



E m is the complex energy eigenvalue, £ m = (h 2 /2fi)k 2 n , and the constant factor Ci{k m ), multi- 
plying <f>i(k m ,r) in equation (|0|), is 



Ciikrn) 



2fi 1 

h 2 2k 



k 



+ 111 df e (k m ) l(d 2 h{-k')^ { d?h{-k') 



2 fe(k m ) dk m 6 V dk' 2 



dk' 3 



(41) 



The normalization constant N^i is now 



2/i 



N m i — \ — 



h 2 ) 16ifc2f+3 



fl{k ri 



d 2 M-k>) 
dk' 2 



(42) 



The expression ( |3"9"| ) suggests the following normalization rule for the chain of Gamow- Jordan 
generalized eigenf unctions belonging to a double zero of the Jost function 



(k mi ' ) 
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4>e(k m ,r) 



(43) 



and 



mi ' J 

Substitution of © and © in © gives 



mt 



-^ e (k m ,r). 



(44) 



2iriRes 



-<^\4 + \k') ><4 + \k')\ x > 

7T 



< $|« m £(fc m ) >< u mi (k m )\x > + < $|u m<? (/c m ) >< u m £(fc m )|x > 
where, the notation means, 



(45) 



<®\ume(k m )>= I <&*(r)u me (k m ,r)dr 



and 



(46) 



< Uml{k m )\x > 



u m e(k m ,r) is defined in (|44|). 

Finally, substitution of the expressions 



{k m ,r)x{r)dr, 



(47) 



and (55) in (^Bp gives the following expansion 



< $|x > = ^ < >< u^lx > + < $I*W >< Unt\X > 

s bound states n^m resonances 

+ < ®\u me (k m ) X u m£ (k m )\x > + < <§>\u m i{k m ) X u ml ik m )\x > 



+ 



7T 



(48) 



This expression shows that, when the Jost function has many simple zeroes and one double zero 
in the fourth quadrant of the complex k-plane, the Gamow eigenfunctions u n i(k m , r) associated 
to simple zeroes of the Jost function and the chain {u m i(k m , r), u m t(k m , r)} of Gamow- Jordan 



generalized eigenfunctions pq, pq , p7[ associated to the double pole of the Jost function are 
basis elements of an expansion in generalized bound and resonant state eigenfunctions plus a 
continuum of scattering functions of complex wave values k'. 

Omitting the arbitrary function $(r) in (f5J), we obtain the complex basis expansion of an 
arbitrary square integrable and well behaved function x(r) 



s bound states n^m 

+ u m£ (k m ,r) < u ml \x > +u m£ (k m ,r) < u ml \x > 
+ - [4 + \k',r) <4 + \k')\ X > dk'. 



(49) 
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In this expression u n i(k n ,r) are the Gamow eigenfunctions representing decaying states 
associated to simple resonance poles of the scattering wave function ipg(k, r), the matrix S(k) 
and the Green's function G^ + '(k;r,r'). The set {w m ^(^m, r), u m i(k m , r)} is a Jordan chain of 
lenght two of generalized Gamow- Jordan eigenfunctions associated to the double pole of the 
scattering matrix S(k) and the Green's function G^(k;r,r') at k = k m . The last term in the 
right hand side of (fS|) and ([i9|) is the background integral defined along the integration contour 
shown in Fig 1. 



7 Jordan blocks in the complex energy basis 

Once it has been established that the Gamow eigenfunctions u n e(k n ,r) and the Jordan chain 
{u m i(k m , r), ii m i(k m , r)} of generalized Gamow- Jordan eigenfunctions are elements of the basis 
set of eigenfunctions in the expansions ([48|) and (|4*9|), we may represent any operator f(Hr^), 
which is a regular function of the Hamiltonian Hr , in terms of its matrix elements in this 
basis. 

Let us start by deriving an expression for the action of f(Hr^) on the generalized Gamow- 
Jordan eigenfunction u m t{k mi r) . With this purpose in mind, let us write the eigenvalue equation 
satisfied by u m e(k m ,r) as 



H^u me (k m , r) = £ m u mE (k m , r), (50) 

where, 



m = - 



h 2 



cP 



2/x Idr 2 



— v[r) — 



(51) 



v(r) is a well behaved short ranged potential which satisfies the conditions stated in section 1. 
Now, let us consider a holomorphic function f(£) of the complex variable £, such that, 



(52) 



3=0 



the coefficients aj are independent of £. 
Then, from (M) and (M), 



f(H^ e) )u me (k m , r) = f(£ m )u mi (k m , r). (53) 
Taking derivatives with respect to the eigenvalue £ m in both sides of (|53|), we obtain, 



j^jj(i)^ u me(k m , r) 



d£ n 



f( c ,du me (k m ,r) df(£ m ) 

J {£rn> ^ 1 ^ U ml {k m , r J 



d£ n 



d£ n 



(54) 
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From this equation and the definition, equations (pEOP , (f|Tj) and (pH]), of u m £(k m ,r), it follows 
inmediatly that 



f(H^)u mi (k m ,r) = f(£ m )u m e(k m ,r) 



d£ 



(55) 



Notice that a necessary and sufficient condition for the existence of du m g(k m ,r)/d£ m is the 
vanishing of ( df (— k)/dk J 

The rule stated in equation (55) permits us to calculate the action of f(Hr) on the gener- 
alized Gamow- Jordan vectors occurring in the complex basis expansions (|48|) and ([491) . 



Now, we can write the operator f(Hr) in terms of its matrix elements in the complex 
energy basis. This may be done by acting with f(Hr^) on the left in both sides of equation 



(ED 



f{Hl l) )x(r) = f(£s)v s i{r) < v sl \x > + ^ f( E n)u n e{k n , r) < u nl \x > 



+ 



(j{£ m )u m t{k m ,r) + u mt >(k m ,r)^ < u mi \x > 



+ f{£m)u m e{k m ,r) < u ml \x > 



+ 



f{E')<j>f\k\r)<^\k')\x>dk'. 



(+)< 



(56) 



Multiplying both sides of Q5TS|) by $*(r) and integrating over r, we get, 
<Q\f(H®)\x> = J2<®\v se > f(£ s ) <v se \ X > <®\u ne > f(£ n ) <u ne \ X > 

s njtm 
+ < <&\u m £ > f(£m) < U mi \x > + < ®\u mi > (f{£ m ) < Urrd\X > 

+ < ume\x >)+-/< *\4 +) W > f{£') < 4 + \k')\x > dk'. 

(57) 

To simplify the notation, suppose that the system has no bound states only resonances and 
that the first two resonances are degenerate. Rearranging equation Q57D in matrix form, we get 



< $\f(H<Q)\x >= (< $\uu >, < $\uu >, < ®\U3i >, 

fm) 



X 



o f(£i) 



o 



V • 
















. 


A 


/ < u u \x >\ 











. 




< u u \x > 


m) 








. 




< u 3e \x > 





f(£*) 





. 
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\ ) 



+- / < mf\k') > f{£>) < ^ +> (k')\ x > dk'. 



(5? 
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In this matrix representation of f(Hr^)f\, the upper left 2x2 submatrix is a Jordan block 
27 1 associated to the chain of Gamow- Jordan generalized eigenf unctions 
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of rank two| 

{uu(ki, r), Uu(ki, r)} belonging to the double zero of the Jost function fi(—k) (double pole of 
the scattering matrix and the Green's function). Except for this 2x2 block, this matrix is 
diagonal with the eigenvalues f(£ n ) in the diagonal entries. Simple zeroes of the Jost function 
correspond to simple (non-repeated) eigenvalues of f(Hr ) while the double zeroe of fe(—k) 
correspond to the twice repeated (degenerate) eigenvalue f{£\) occuring in the Jordan block. 
The off-diagonal non- vanishing element in this block is df{£i)/d£\. 

The difference in physical dimensions of the off-diagonal and the diagonal entries in the 
2x2 Jordan block is compensated by the difference in normalization of the Gamow- Jordan 
chain {uu(ki,r),uu(ki,r)} and the Gamow eigenfunctions u n i(k n ,r) (n = 3,4,....) which are 
normalized according to fl4"3] , and ( |T7| , [H| |3"3| ) respectively. 
It will be instructive to consider some simple examples. 



We first choose f{H, 



Hr . Then, from 



we obtain, 



< 



X 



x 



(Ex 


1 

















Si 




















£3 




















£4 









/ < u u \x >\ 

< uu\x > 

< Uai\x > 



V / V • / 

+- [< $\4 + \ k ') >E'< 4 + \ k ')\x > dk'. 

7T Jc 



(59) 



From this example, it is evident that in a degeneracy of two resonances in the absence of 
symmetry, the degenerate complex eigenvalue £\ occurs twice in the spectral representation 
of the radial Hamiltonian Hr given in (|59D, while there is only one Gamow eigenvector or 
normal mode, uu(ki,r), associated to the degeneracy. This is so, because the Gamow- Jordan 
generalized eigenfunction or abnormal mode, uu(ki,r), is not an eigenfunction of the radial 
Hamiltonian Hr . This is a generic property of this kind of degeneracy which may be stated 
in slightly more formal terms as follows: In a degeneracy of resonances in the absence of 
symmetry, the algebraic multiplicity is always larger than the geometric multiplicity. Here, we 
mean by algebraic multiplicity of a degeneracy, /i a , the number of times the degenerate complex 
eigenvalue is repeated, and, by geometric multiplicity of the degeneracy, fi g , the dimensionality 
of the subspace spanned by the eigenvectors associated to the degenerate eigenvalue p5|, |26| , ^7j| . 

Then, 



Ha > fl g . (60) 

""^From the way it was derived, it is evident that the matrix in equation (^8|) represents the action of f{Hr ) 
as an operator on the space of continuous antilinear functional on the Schwarz space of very well behaved test 
functions. 
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Let us consider now, the complex energy representation of the resolvent operator. In this 
case f(Hr ) = 1/(E — Hr ). Then, from (|5lf ), we obtain, 



< $1 



1 



E-H, 



(0 



\x 



X 



/ 1 1 

' E-Ei (E-Etf 

1 








E-Ex 










E-£ 3 











E-£ 4 















/ < u u \x >\ 




< uu\x > 




< u 3£ \x > 


• •/ 


\ J 



j,<4 + \k')\ X >dk'. 



(61) 



It may easily be verified that, when we delete the arbitrary functions $(r) and x{ r ) m 
this expression, the resulting expansion for < r\ — 1 w \r' > is just the expansion in resonance 
eigenf unctions of the complete Green's function 



G ( e + \k;r, r') 



+ 
+ 



2/2 



E 



v si (k,r)v* e (k,r') 



bound states 



E+\E S \ 
u n e(k n ,r)u ne (k n ,r') 



E — E „ 



E 

nj^m resonant states 

r)u m e{k 

mi 1 I 

{k m , r)u m e(k m , r') + u m e(k m , r)u m e(k m , r') 



2 

+ - 



7T 



c 



(E — £ m ) 
(k 2 - k" 2 ) 



dk'. 



(62) 



The occurrence of the double pole in G[ + \k; r, r'), as function of the complex energy, is thus 
associated to the occurrence of a Jordan block of rank two in the complex basis representation 
of the resolvent operator and a Jordan chain of Gamow- Jordan generalized eigenfunctions 
{uu(ki,r),uu(ki,r)} associated to the double zero of the Jost function. 

Finally, let us consider the time evolution operator exp (—iHt). For each fixed value of the 
angular momentum, it will be enough to consider the operator f(Hr^) = exp (—iHrt). In this 
case, from equation (EH) 
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f exp (— i£it) 




V 



< $| exp (-iH^>t)\x >= (< ®\u u >, < $\uu >, < $\v>M > 

itexp(-iSit) . A /< uu\x >\ 

exp (—i£it) . . < uu\x > 

exp {-i£ 3 t) 

exp (—iS^t) 



< u 3i \x > 



I \ 



J 



+- / < > exp (-i£'t) < 4 + \ k ')\x > dk'. (63) 

7T 

As in the previous examples, the time evolution operator is non diagonal in the complex 
energy basis representation. The time evolution of the Jordan chain of Gamow- Jordan general- 
ized eigenfunctions {uu(ki, r),uu(ki, r)} is given by a Jordan block of 2 x 2 with an exponential 
time dependence in the diagonal entries and a first order polynomial times an exponential in the 
off-diagonal entry. Hence, the time evolution of the Gamow- Jordan generalized eigenfunction 
or abnormal mode is a superposition of the abnormal mode uu(ki,r) evolving exponentially 
in time plus the normal mode Uu(k,r) evolving according to the product of a first order poly- 
nomial times an exponential time evolution factor. The time evolution of the normal mode 
uu(k\,r) in the Gamow- Jordan chain {uu(ki, r), uu(ki, r)}, as well as the time evolution of 
all other normal modes u n e(k n , r) associated to the simple zeroes of the Jost function (simple 
poles of the scattering matrix) is purely exponential. 



8 Orthogonality and normalization integrals for Gamow- 
Jordan eigenfunction. 

As in the case of bound and resonant state eigenfunctions associated with simple poles of the 
Green's function, we may derive orthogonality and normalization rules for the Gamow- Jordan 
eigenstates in terms of regularized integrals of the generalized Gamow- Jordan eigenfunctions. 



Following the same procedure as in Berggren|23|, pMj , it may be shown that, when fe(-k') has 
a double zero at k! = k m , the following relations are valid, 



and 



1 



i8k m 



d 2 f e (~k>) 
dk' 2 



lim 



r 2 d(j)j(k m , r) 
dk m 



(f> e (k m ,r)dr 



(64) 



^Sk^ +1) 



fe{K 



l( d 3 fe(-k') 
L3V dk' 3 



d 2 f,{-k') 
dk' 2 



lim 



dfe{k r , 



fi(k m ) dk, 
d(j) e (k m ,r^^ 2 



dr. (65) 



From the expression (|41|) for Ce(k m ) and equations (0) and (|65"D, it follows that, 
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2Q(*».) — [~ e- ^'^' r) MK>,r)dr) = 0, (66) 
u V-»o /„ dk m J 



/i V->o 

which may be rewritten as, 



u,ft, m 



lim 

v-±0 



d<pe{k rt 



d£r, 



+ C e (k m )(j>e(k m , r) dr = Cj(k m ) lim 



v->0 



>f(fc m ,r)dr, (67) 



but, according to equation fl22j ) and (p3|) , when ft(—k) has a double zero at = fc m , the integral 
in the right hand side of (|67|) vanishes. Therefore, the integrand in the left hand side of ( |67|) 
is the square of the generalized Jordan- Gamow eigenf unction and the relation ( |65D translates 
into 



oo 

vr 2 12 1 



lim / e~ ur (j)f(k m ,r)dr = (68) 







which shows that also the regularized integral of the square of the generalized Gamow- Jordan 
eigenfunction vanishes. 

An expression for the normalization constant A/^ in terms of a normalization integral may 
be obtained from (§54]), 



2d4> e (k m} r) 



A/^ = lim/ e— ^ tK ': m, ' J Mkm,r)dr, (69) 



writing d(pi/d£ m in terms of (j>e(k m ,r) and recalling that the integral of <fij(k m ,r) vanishes, we 
get, 



M 2 ml = lim / e~ ur 2 Mk m , r)Mk m , r)dr, (70) 
v -° Jo 

which shows that the right hand side of ( [70]) is the normalization integral for the Gamow- Jordan 
generalized eigenf unctions associated with a double pole degeneracy of resonances with A/j^ as 
given in fl42|) . However, it is convenient to note that this expression does not fix the normaliza- 
tion rule for (j>i(k m ,r) and (f>t(k m ,r) in a unique way. Since <f>i(k m ,r) and (j>e(k m ,r) are linearly 
independent, they have different dimensions and its product has no obvious interpretation in 
terms of observable quantities, therefore, there is no a priori reason to normalize both functions 
with the same normalization constant. Thus, we still have the freedom to write (|70|) as 



lim f°° e— 2 (^Mk m , rj) (—l—^km, rj)dr = l, (71) 
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where A/^ is given in fl4"2"|) and X m is a non-vanishing real or complex number that we asso- 
ciate with the double pole singularity of Gg(k; r, r') at k — k m . Therefore, a more general 
normalization rule for the Gamow and Gamow- Jordan generalized eigenfunction that the one 
proposed in (f4|), © and fl4|) would be 



mi ' I 



-4>e(k r , 



(72) 



and 



mi ' I 



X 

-^(f) e (k m ,r) 



(73) 



' ml 



With this normalization, the orthogonality and normalization integrals for the generalized 
Gamow- Jordan eigenfunction associated to a double pole of the Green's function, equations 
(B8|) and (|70D, take the form 



lim 



'u 2 ml (k m ,r)dr = 



(74) 



lim 



'u 2 mt {k mi r)dr = 



(75) 



and 



lim / e vr u me (k m ,r)u mi (k, r)dr = 1. (76) 
Jo 

The form of these orthogonality and normalization conditions is independent of the value 
of the constat X m . However, if the Gamow- Jordan generalized eigenfunction are normalized 
according to (^) and flT3"|) the expression for the residue at the double pole of Gg (k;r, r') 
would be explicitly dependent on X m , since a factor X 2 m will appear multiplying the term 
u m e(kmi r )u m e(k m , r') in the expression for the residue at the double pole of G^\k; r, r') given 
in equation flB"2"|). 



-^-m u rne(kmi r )um£(k rnj r') u m f>(k m , r)u m f (k m , r') + u m j>(k m ,r)u m g(k m} r') (77\ 

(e - s m y + (e - £ m ) • 1 ) 

As is evident from the definition (0), the generalized eigenfunctions (p n e(k m ,r) and 4> n £(k m ,r) 
have different dimensions, if one takes X m of dimension (energy) 1 / 2 the normalized eigenfunc- 
tions u n i(k n ,r) and u n i(k ni r) have the same dimensions namely (energy)" 1 / 2 so that when 
(X m ) = (energy) 1 ^ 2 the higher order Gamow- Jordan vectors become Jordan vectors with the 
same dimensions as the Gamow vectors. 

This freedom in the normalization rules could be used to define normalized Gamow- Jordan 
eigenfunctions with the same dimensions as those of the Gamow eigenfuntions associated to 
simple poles of G^\k; r, r'). 
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9 Summary and conclusions 



In the theory of the scattering of a beam of particles by a short ranged potential, resonances 
are asociated to the occurrence of poles of the scattering matrix Se(k), the Green's function 
(k;r, r') and the scattering wave function if>ne(k,r). These resonance poles are caused by 
zeroes of the Jost function lying in the fourth quadrant of the complex k— plane. Accordingly, a 
degeneracy of resonances, that is, the exact coincidence of two (or more) simple resonance poles 
of the scattering matrix, results from the exact coincidence of two (or more) simple resonance 
zeroes of the Jost function, which merge into one double (or higher rank) zero lying in the 
fourth quadrant of the complex A;— plane. 

We found that, associated to a double resonance zero of the Jost function, there is a Jordan 
chain of length two|25|, |26|, |27|] of generalized Gamow- Jordan eigenf unctions {?w(/c m ,r), 
u m i(k m ,r)} belonging to the same degenerate complex energy eigenvalue £ m . Hence, the cor- 
responding second rank pole occurring in the scattering matrix, Se(k), the Green's function, 
G^~\k; r, r'), and the scattering wave function, ip^ + \k, r), is also associated to this Jordan chain 
of Gamow- Jordan generalized resonance eigenfunctions. 

As the two simple zeroes of the Jost function merge into one double zero, the two Gamow 
eigenfunctions corresponding to the two resonances that become degenerate merge into one 
Gamow eigenfunction or normal mode belonging to the double zero of the Jost function. The 
other element in the Jordan chain, namely, the Gamow- Jordan generalized eigenfunction or ab- 
normal mode is not a proper eigenfunction of the radial Hamiltonian. Hence, at a degeneracy 
of resonances, one resonance eigenfunction or normal mode is lost, and a new kind of gener- 
alized resonance eigenfunction or abnormal mode is generated. Therefore, the dimensionality 
of the subspace of eigenfunctions associated to a degeneracy of two resonances or geometric 
multiplicity, fj, g , of the degeneracy is one, yet, the number of times the degenerate complex 
energy eigenvalue is repeated in the spectral representation of Hr or algebraic multiplicity of 
the degeneracy, fi a , is two. It follows that, the algebraic multiplicity is larger than the geometric 
multiplicity of a degeneracy of resonances. 

Explicit expressions for the normalized Gamow and Gamow- Jordan generalized eigenfunc- 
tions in the Jordan chain, written in terms of the outgoing wave Jost solution, the Jost function 
and its derivatives evaluated at the double zero, are obtained from the computation of the 
residue of the scattering wave ip^ + \k, r) function at the double pole. 

We also showed that the Jordan chain of generalized eigenfunctions are elements of the 
complex biorthonormal basis formed by the real (bound states) and complex (resonance states) 
energy eigenfunctions which can be completed by means of a continuum of scattering wave 
functions of complex wave number. With the help of this result, we derived expansion theo- 
rems (spectral representations) for operators f(Hr) which are regular functions of the radial 
Hamiltonian Hr^ ■ In this basis, the operator f(Hr^) is represented by a complex matrix which 
is diagonal except for one Jordan block of rank two [p5| , |26| , [27j associated to the double zero of 
the Jost function and the corresponding chain of generalized eigenvectors. The diagonal entries 
in this matrix are the eigenvalues f(£ n ), simple zeroes of the Jost function correspond to non- 
degenerate eigenvalues of f(Hr ) while the double zero of the Jost function corresponds to the 
twice repeated (degenerate) eigenvalue f(£ m ) in the diagonal entries of the Jordan block. The 
off-diagonal, non-vanishing element in this block is df(£ m )/d£ n . In particular, the occurrence 
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of a double pole in the Green's function, as function of the complex energy, is thus associated to 
the occurrence of a Jordan block of rank two in the complex basis representation of the resolvent 
operator and the corresponding Jordan chain of Gamow- Jordan generalized eigenfunctions. 
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